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Abstract. - We report experimental and numerical results on the growth of a super stable 
heap (SSH). Such a regime appears for flows in a thin channel and for high flow rate : the flow 
occurs atop a nearly static heap whose angle is stabilized by the flowing layer at its top and 
the side wall friction. The growth of the static heap is investigated in this paper. A theoretical 
analysis inspired by the BRCE formalism predicts the evolution of the growth process, which 
is confirmed by both experiments and numerical simulations. The model allows us to link 
the characteristic time of the growth to the exchange rate between the "moving" and "static" 
grains. We show that this rate is proportional to the height of the flowing layer even for thick 
flows. The study of upstream traveling waves sheds new light on the BCRE model. 



Introduction. - The dynamics of granular media appears to be very different from that 
of usual fluids [1]. In particular such media have the ability to remain static (i.e. a solid- like 
behavior) at a non-zero tilt angle. When this tilt angle exceeds a critical value a surface flow 
is observed : most of the pile remains static except for a relatively thin rolling layer [2,3]. A 
granular medium can thus display both solid-like and fluid-like behaviors. Very recently [4], it 
has been shown that when the granular medium is confined between two vertical glass walls, 
the angle of the free surface increases with increasing input flow rate. Unusually steep granular 
heaps (Super Stable Heap - SSH) can then be obtained. This letter presents an experimental 
and numerical study on the growth of such heaps which allows one to characterize the exchange 
rate between the fluid-like and the solid-like phases. 

Summary of the previous work on Super Stable Heaps. - All experiments presented in 
this paper were performed in the same channel as in our previous work [4] . Our experimental 
set-up consists of a thin three-dimensional open channel that could be inclined from horizontal 
to an angle 9 (figure ^ shows a similar system used in our simulation) . Although 9 could 
range from 0° to 90°, all experiments were conducted in a steep configuration in which 9 was 
always set to be higher than the angle at which all the grains flow out of the channel after 
the feeding was cut off. Experiments were performed using Fontainbleau sand (round shaped, 
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0.3-0.5 /im). Granular material was continuously poured on a rough bottom placed between 
two transparent vertical glass plates through a hopper whose aperture precisely controlled the 
input flow rate Qi n (mass per unit of time and per unit of width). The rough bottom was 
made of grains glued to a flat rigid substrate. The gap between the two vertical plates W was 
varied from 5 mm up to 20 mm and the length of the channel, L, could by varied by moving 
the hopper. For further details see [4]. 



Fig. 1 - a) (In color online) Snapshot of the flow from numerical simulation. Velocity color-scale b) 
Sketch of the transient regime. The dark part corresponds to moving grains, the gray part to static 
grains. 

For low flow rates we observed classical flows on the rough bottom of the channel but a 
new regime, occuring for high flow rates, was investigated. In this regime, a static wedge-like 
heap gradually forms: some of the flowing grains get trapped and contribute to the growth 
of the static heap. Let us mention here that the grains are not completely static but can 
experience a slow creeping motion. For simplicity, these grains will however be refered to 
as static. The system eventually reaches a steady state consisting of a uniform flow atop a 
static heap. Let us mention here that although the flow is uniform when averaged over time, 
there can exists local fluctuations of the flow height, h, sometimes corresponding to upstream 
traffic waves described in [5] . Following the definition in [4] , h was measured using the velocity 
profiles. When the feeding is cut off, the heap slowly erodes and falls apart. This shows that 
the heap is dynamically stabilized by the flow at its top. One remarkable feature exhibited by 
such heaps is their unusual steepness (up to 60° in our experiments) and this is why they are 
referred to as Super Stable Heaps. Figure ^ is a sketch of the regime described here. The 
SSH appears in light gray and the flowing layer in black. For given values of Qj„ and W, the 
system "chooses" an angle ip for the free surface and the corresponding height of the flow h. 
The values of ip and h depend only on the input flow rate, Qi„, and on the channel width, 
W, and not on the tilt angle 9 nor on the channel length, L. A simple depth-averaged model 
provided us with a simple law linking the angle of the SSH in the final steady state, the height 
of the flow and the channel width. Writing the momentum balance leads to: 




h 



(1) 



tan ipoe = [i t 




where fa and fi w are two positive material constants. 
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Our previous study focused on the steady state and identified the mechanisms responsible 
for the existence of SSH. The present paper focuses on the growth of the SSHs i.e. the time 
transient leading to the steady state. The study of this transient regime is of great interest on 
its own and can moreover provide one with quantitative information on the exchanges existing 
between static and moving grains when granular material is flowing on a static heap. Our 
experiment might then shed some light on the BCRE model [6] . 

Preliminary Observations. Before discussing the results let us mention two simple 
observations that simplified the detailed study of the growth of the heap. First, we observed 
that the growth was a "slow" process: among the many grains entering the channel, only very 
few of them actually contributed to building the heap. Figure|5|is a plot of m out , the mass of 
grains having exited the channel, vs. time and shows that the output rate Q ou t slowly tends 
toward a constant. This is clear on the inset which is a plot of Qout — dtm out . One can 
see that even after 30 seconds, a period much greater than the time it takes a grain to flow 
through the channel (typically Is), the steady regime is not yet reached. Let us mention that 
the variation of the height, h, was found to be very small during the growing process. 
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Fig. 2 - Plot of the output mass m out vs. time measured experimentally. Inset: Qout = dtm vs. time 



The second preliminary observation is that the surface of the heap remains flat throughout 
the whole duration of the growth. That is to say that the heap maintains a linear profile while 
its slope changes over time. One can therefore define a single angle for the whole heap, ip t , at 
any given time t (see figure H 3 ). 

The free surface presents no curvature in the z direction, meaning that a side-view is 
sufficient to describe the system. For technical reasons the surface area of the heap, A t , was 
easier to measure than the angle tp t . In both our experimental system and our simulation, 
the channel was long enough to allow for the neglect of the surface area of the flowing layer 
relative to that of the static heap (less than 5%). The surface area of the whole heap (SSH 
and flowing layer) was then assumed to be equal to the surface area of the SSH. This allowed 
for easier measurements of A t . Given these observations, a simplified geometrical expression 
for the surface area of the heap can be given, as long as iff — 9 is not too large (i.e. j 30): 

Theory. - As mentioned above, the flow consists of a layer moving rapidly atop a static 
heap. Since a flowing grain can be trapped and conversely a static grain can start moving 
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after a collision, there exist exchanges between the two species of grains. The distinction 
between two species of grains is the basis of the BCRE model [6]. In this model, the height 
profile of the flowing layer obeys a convective diffusion equation, and the exchanges between 
the two species are ruled by a local law: accretion is proportional to the difference between 
the neutral angle (i.e. the angle at which erosion compensates accretion) and the local slope 
of the flow. The positive proportionality constant, denoted T is called the exchange rate and 
is proportional to the flow height, h [6]. However, it was suggested in the literature [7] that 
this should not hold for thick flows (h > 5-10 d), for which T should be proportional to v up , 
the velocity of upstream traveling waves, the existence of which has been reported in the 
literature [5]. 

Our idea was to use a law similar to that of the BCRE model in order to describe the 
overall growth rate of the SSH. In our experiments, the neutral angle was none other than 
tfoo, the angle of the SSH in the final steady state since then the angle of the heap remains 
constant, meaning that erosion exactly compensates accretion. Therefore, we assumed that 
the surface area of the SSH follows the law: 

d t A t = v{ Voa - tpt) , (r > o) (3) 

Using equation and leads to: 

L 2 

A t = (1 - e _ r), where r = — (4) 

This model predicts an " 1-exp" growth and provides us with a characteristic growth time, 
which appears to be inversely proportional to the exchange rate, T. 

Simulation methods. - Numerical simulations were performed using a Molecular Dynam- 
ics (MD) method, a.k.a. Discrete Element Method (DEM). This type of simulation has been 
widely used in the past and has proven to be a very reliable method [8]. Figure^, presents 
a snapshot of the flow obtained with our simulation. In the simulation, we used a dashpot- 
spring force model : F[j — k n Sij + "f n dtSij for the normal force and a Coulomb friction law for 
the tangential force : — min(fj,F^ , 7tV^-) , where 5^ is the overlap between two colliding 
particles, k n a spring constant, j n a viscous damping constant, fi a friction coefficient, 
the sliding velocity of the contact, j t a visous regularization constant. The following values 
were used : particle radius R—A mm, mass=0.16 mg, k n = 40000 N.m , j n = 1.2 s _1 and 
7t = 5 s" 1 , leading to a coefficient of restitution of 0.4. Different values of the latter were tried 
but did not affect the flow properties, which was expected since in dense assemblies of grains 
the effective coefficient of restitution is close to zero, regardless of the material properties. 

Because of the nature of the phenomenon, it was not possible to use periodic boundary 
conditions. The channel is constantly fed by releasing grains from a virtual hopper. The 
rough bottom of the channel consisted of immobile particles. The collisions against the side 
walls were treated like particle-particle collisions (identical simulation parameters n, k n , j n 
etc) with one of the particle having infinite mass and radius, which mimics a large flat sur- 
face. The granular material was made slightly polydisperse (±20% in size) in order to avoid 
crystallization. Our simulations typically contained a very large number of particles, between 
10 5 and 10 6 , depending on the size of the channel and on the input flow rate Qi n and were 
run for typically 10 7 time steps. 

Experimental and Numerical Results. The time evolution of the surface area, A t , 
was studied both experimentally and numerically. In the experiments, A t was measured 
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using imaging software. As for the simulations, A t was computed from the number of grains 
belonging to the SSH. 




Fig. 3 - Surface area At vs. time. Insets: semi-log plot of (A^ — At)/ Aoo- a): experiments, solid line: 
L = lm, W = 5mm, Q — 2.4 gr.s -1 . mm -1 , dashed: L = lm, W = 10mm, Q = 11.1 g. s _1 . mm -1 . 
b) numerics, solid line: L = 200 R, W = 10 R, Q = 225 part. s" 1 .mm" 1 , dashed: L = 200 R, W = 
20 R, Q = 650 part. s' 1 .mm' 1 

Figure shows plots of the surface area of the SSH vs. time measured experimentally and 
numerically for various values of the control parameters L, W and Qi n - The evolution seems 
to follow an 1-exp law, which is confirmed by the right-hand side of the figure. The insets on 
Figure show semi-log plots of {A^ — At) /Aoo vs. t and presents excellent agreement with 
equation since all data follow a linear law. Indeed, both experimental and numerical data, 
for a wide variety of parameters, present a remarkable agreement with the "1-exp" law. One 
can therefore measure r and compute the value of the exchange rate, T, for various values of 
the control parameters. 

In this last section, we would like to comment on the exchange rate T. As mentioned above 
the value of T can be extracted from the characteristic time, r, using: r = L 2 /(2T). The 
L 2 scaling was expected since the properties that were monitored in this study are surface 
areaand checked numerically. A more relevant, geometry-independent exchange rate, 7, can 
then be defined: 7 = 2T/L 2 = 1/r. Measurements of r were performed for various values of 
W and Qi n , allowing for the study of r as a function of the angle (p^. Figure^ is a plot of 
r vs tan ifoa. The straight lines do not come from a fit but only help to visualize the general 
trend. 

The first conclusion that can be drawn from this plot is that r depends on the channel 
width W . That is to say that for a given angle <poo, r varies when W is changed. Namely, 
the thinner the channel, the larger the characteristic time. In our previous work on SSHs we 
showed that there exists a simple relationship between ^oo, h and W (equation QJ). In the 
light of this previous result, one can interpret the data in a new fashion. For a given ip 00, T 
decreases with increasing W, as for a given W, r decreases with increasing <poo, meaning in 
both cases that 7 increases with increasing h. The function t(</? 00 , W) can be expressed as the 
new function "f(h), which is easier to interpret. As mentioned above, the heigth of the flow 
remains constant during the growth, and can therefore be computed as follows. The values 
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Fig. 4 - a) Characteristic growth time r vs. tan tpoo measured experimentally, b) exchange rate 7 vs. 
flow height, h 



of the coefficients \ii and /i w having been determined in our previous work (/i^=tan 23.3° and 
/i„,=tan 33.7°), one can then compute the flow height in the steady regime of each run using 
equation £3>. One can then plot 7 as a function of h. Figure^]} presents such a plot and shows 
a decent data collapse. 7 is found to be a linear function of the flow height h. In agreement 
with theoretical work based on the BCRE model [3], the exchange rate is proportional to the 
height of the flow but in our case, no saturation was observed. 

Although 7 is clearly a linear function of h, it remains interesting to check whether 7 
can be expressed as a function of v up , the velocity of upstream traveling waves. In order to 
compute this velocity, we measured the fluctuations of the height, h(x,t), in the final steady 
state. Figure is a grey-scale space-time diagram of h(x,t) — h. On this figure, a dark spot 
indicates a local height, h(x,t), greater than the average i.e. a hump, whereas a white spot 
corresponds to a depression. 




Fig. 5 - Space-time diagram of the fluctuations of the height of the flowing layer, h(x,t) — h. 

Two patterns can be observed on figure [SJ: on one hand, quasi- vertical narrow lines, on the 
other hand, wide oblique stripes. The former correspond to perturbations "floating" at the 
free surface (like a object floating on a river) and therefore moving downward at the velocity 
of the free surface, v sur f aC e, whereas the latter corresponds to upward traveling waves. The 
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main direction of both kinds of stripes (i.e. v sur f ace and v up ) can be extracted using an 
angular FFT method or directly measured by fitting them with straight lines. The narrow 
lines (corresponding to downward motion) appear to be less oblique than the wide stripes 
(corresponding to upward waves), meaning that v sur face > v up (on figure El : v surface— 
82 cm/s and v up =11.8 cm/s). In other words, the upward travelling waves travel slowly 
compared to the characteristic velocity of the flow. Secondly, the spacing between the wide 
stripes seems to be regular, indicating that there exists a well-defined wave-length. The origin 
of this wave-length is unclear and might be related to the phenomena responsible for the 
onset of the upward travelling wave. The velocity v up was measured for different values of the 
control parameters but no simple relationship could be found between r and v up . This last 
result seems to indicate that even for the thick flows we observed (h ~ 50d), the exchange 
rate is not proportional to v up . Let us mention here that the use of v up in the BCRE model 
has already been questioned in the literature [9] . 

Conclusion. - The growth of the Super Stable Heap was studied both experimentally 
and numerically (using the Molecular Dynamics method) . Among the many grains flowing at 
the surface of the static heap, a few get trapped and contribute to the formation of the heap. 
During the growing process the free surface of the flow was found to remain flat. This allowed 
us to introduce a model for the accretion/erosion process, similar to the BCRE model [6]. 
A global accretion/erosion law predicted the time evolution of the surface area of the heap 
that was confirmed by both experiments and simulations. The model provides one with a 
characteristic time which is a direct measurement of the exchange rate, 7. We showed that 
7 is proportional to the height of the flow, h, even for thick regimes (h ~ 50<i). However, no 
simple relationship could be established between 7 and v up , the velocity of upward traveling 
waves (traffic waves). It was suggested in [7] that they should be proportional for thick flows. 
This disagreement might indicate that 7 and v up are never proportional or that there exists a 
fundamental difference between flows on SSHs and on classic heaps. The relationship between 
7 and v up should be investigated experimentally in the latter case. 
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